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Abstract 

The problem of constructing internally rotating solitons of fixed angular frequency 
ui in the Faddeev-Skyrme model is reformulated as a variational problem for an energy- 
like functional, called pseudoenergy, which depends parametrically on uo. This problem 
is solved numerically using a gradient descent method, without imposing any spatial 
symmetries on the solitons, and the dependence of the solitons' energy on u), and on 
their conserved total isospin J, studied. It is found that, generically, the shape of a 
soliton is independent of u), and that its size grows monotonically with u. A simple 
elastic rod model of time-dependent hopfions is developed which, despite having only 
one free parameter, accounts well for most of the numerical results. 

PACS classification numbers: 05.45.Yv, 03.50.-z 
Keywords: topological solitons, knot solitons 



1 Introduction 



Many field theories of interest in fundamental physics support topological solitons - spatially 
localized, stable lumps of energy whose strongly particle-like characteristics make them nat- 
ural theoretical models of elementary particles. Perhaps the best developed model from this 
viewpoint is the Skyrme model, whose solitons are posited to model atomic nuclei. It is of 
fundamental importance in this context that individual solitons possess both rotational and 
internal rotational (or isorotationat) degrees of freedom. In the Skyrme model, the rota- 
tional degrees of freedom account, after quantization, for the spin of atomic nuclei, while the 
isorotational degrees of freedom account, roughly speaking, for their difference in "flavour". 
For example, both protons and neutrons are modelled by degree one solitons, the only dif- 
ference being in the sense in which they are (internally) spinning. In practice, finding even 
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static classical solitons is a significant computational challenge, so it is not surprising that 
almost all studies of (iso) spinning solitons work within a rigid body type approximation: the 
(iso) spinning soliton is assumed to have precisely the same shape as the static soliton, but with 
its spatial and internal orientation allowed to vary in time [2, 10]. It has long been recognized 
that this is not a very satisfactory approximation, and various attempts have been made to 
improve on it [14, 11]. With a few exceptions [1, 3], these have fallen short of numerically 
solving the full field equations to construct genuine (iso) spinning solutions. 

In this paper we construct genuine isospinning (i.e. internally rotating) soliton solutions 
in the Faddeev-Skyrme model [4]. Such solitons are conventionally called hopfions, since 
they are classified topologically by their Hopf degree (an integer- valued topological invariant). 
Explicitly, our field if is S 2 valued, lives in M 3+1 , and has Lagrangian density 

c = ■ &<p - \{d^ x d u y) ■ (&><p x er<p) - - cpl). (i.i) 

This model is ideal for our purposes because it is three-dimensional and has a rather rich 
spectrum of static soliton solutions, which should be understood as (possibly linked or self- 
knotted) string-like objects, but has only one internal rotational degree of freedom, because 
the potential breaks the internal rotational symmetry group to 5*0(2). (Of course, one needs 
a potential term in £ in order to allow any isospinning solutions at all.) Hence one can ask 
unambiguously "what is the degree Q isospinning hopfion of given angular frequency u G R?" 
For a given degree Q, there are usually several different stable static soliton solutions of rather 
similar energy (the number of solutions seems to grow with Q; see [5] for a study of static 
hopfions in a model with potential). It is interesting to discover how the energy of such solitons 
varies as their angular frequency (or their conserved isospin) changes. We will find several 
examples where these energy curves cross, so that the most energetically favourable shape of 
the hopfion for a given degree Q changes when one isospins them fast enough. 

Our method is to use the Principle of Symmetric Criticality to reduce the problem of find- 
ing isospinning solitons of angular frequency u to a static "energy" minimization problem, 
where the "energy" to be minimized (more properly called pseudo-energy) F u depends para- 
metrically on u. We solve this problem numerically by a standard gradient descent method. 
Starting at u = (where the solitons are already well understood), we gradually increment 
to and reminimize, in this way constructing curves of isospinning solitons, parametrized by 
u. This method allows us to construct isospinning solitons without imposing axial symmetry, 
in contrast to previous studies which have used simulated annealing [3] or direct solution of 
the field equations [1]. Since very few static Hopf solitons are axially symmetric, this is a 
significant advantage. 

We expect that at least one minimizer of exists for each Q whilever u < min{a;i, U2}, 
where uj\ is determined by the radius of the target two-sphere (here taken to be 1) and U2 = /i 
is determined by the "meson" mass of the model. In fact, one could prove this, at least for 
Q in an infinite (but unknown) subset of Z, by an obvious modification of the methods of 
Lin and Yang [12], who considered the functional Fq with fi = 0. It is well understood that 
solitons cannot spin with u > \i because they become unstable to radiation of mesons. We will 
see evidence, however, that if > 1, solitons generically lose stability and collapse before uj 
reaches /i, due to the extra nonlinear velocity dependence generated by the Skyrme term. This 
is not specific to the model considered here, and should be generic for Skyrme-type models. 
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We are able to get significant analytic insight into the behaviour of isospinning hopfions by 
developing a time-dependent extension of the elastic rod model of hopfions introduced recently 
by two of us, in collaboration with Sutcliffe [7]. This extended model, which is of consider- 
able geometric interest in its own right, predicts a universal energy-frequency relationship, in 
rather good agreement with our main numerical results. We also test the rod model by com- 
paring its predictions for infinite isospinning Hopf strings (i.e. solutions which are equivariant 
with respect to translations in a fixed direction) with numerics. Again, we find rather good 
agreement. 

The rest of this paper is structured as follows. In section 2 we reduce the problem of finding 
isospinning hopfions to an energy minimization problem. In section 3 we develop an elastic 
rod model for time-dependent Faddeev-Skyrme fields and extract phenomenolgical predictions 
on isospinning hopfions and Hopf strings. In section 4 we present our numerical results, which 
are briefly summarized in section 5. 



2 Reduction to a static variational problem 

Since isorotation involves only rotational symmetry of the target space S 2 , it is convenient 
to consider the Faddeev-Skyrme model on a general oriented Riemannian manifold M. This 
allows one to treat in unified fashion the case of principal interest, M = R 3 , and the cases of 
soliton chains or strings, M = R 2 x S 1 , sheets M = R x T 2 , or geometrically nontrivial domains 
(of potential interest for cosmological applications, for example). Given a time- dependent field 
tp : R x M — > S 2 , we have at each fixed time t a mapping tp(t, •) : M — > S 2 which we shall, in 
a slight abuse of notation, again denote tp, and a time derivative tp, which is a section of the 
bundle tp~ 1 TS 2 over M. Using these, we define, at time t, the kinetic and potential energy 
functionals to be 

t = J M \m 2 + \\^(^)\ 2 (2.i) 

v = [ iid^ + Vni 2 + */(¥>), (2.2) 

where Q is the area form on S 2 , ip*Q its pullback to M, i denotes interior product, and 
U : S 2 — > [0, oo) is a smooth potential function which we assume attains its minimum value 
at some point tpoo G S 2 , and is invariant under rotations about ip^. If M is noncompact, we 
assume that tp(t, x) — > tpoo, as x — > dooM, sufficiently fast for all integrals to converge. 

Let oj > be a fixed constant. We seek time-periodic solutions of period 2tt/u. By 
definition, these are critical points tp : x M — > S 2 of the action functional 

S(<p) = [ (T-V) (2.3) 

where = R/(27r/o;)Z is the circle of length 2tt/u. Denote by X w the completion in C l of 
the set of smooth maps tp : x M — > S 2 of finite action. We define an action of the group 
S 1 = R/27rZ on as follows: 

{[a),tp) tp {a] , tp [a] (t,x) = R(a)tp(t- a/u,x) (2.4) 
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where R(a) denotes the SO (3) matrix generating rotation through angle a about the axis ipoo. 
Clearly, S((p[ a ]) = S(ip) for all ([a], <p), since the action is separately invariant under both time 
translation and isorotation about ip^. Denote by X^ 1 the set of fixed points of this action. 
Then ip £ X^ if and only if 

<p(t,x) = R(ut)ip(x) (2.5) 

for some map ip : M — > S 2 . We may think of ip as the stationary field (p when viewed in an 
internally corotating frame. Since S 1 is compact, it follows from the Principle of Symmetric 
Criticality [13] that ip £ X^ is a critical point of S : X^ — > R if and only if it is a critical 
point of the restricted action S : X^ — > R. Now 

S(R(ut)iP(x)) = ^ S^u'Jj^ x + |d(^oo ■ n 2 ) ~ W)} , (2.6) 

so a uniformly isorotating field of the form (2.5) is a critical point of S if and only if the static 
field ip : M — > S 2 is a critical point of the functional 

F u {1>) = J {^(|d^| 2 - a; 2 |d(^oo ■ n 2 ) + l\m\ 2 + (UW - Iuj 2 ^ x ^| 2 )| . (2.7) 

We shall call this functional the pseudoenergy of ip. It has a rather interesting and natural 
form, as we now describe. 

The first two terms of (2.7), taken together, can be interpreted as the Dirichlet energy of 
the map ip : M — )• S 2 , where 5* 2 is given the deformed metric 

(X, Y) u = X ■ Y - u 2 ^ ■ X)(^oo ■ Y) (2.8) 

for all X, Y £ T^S 2 . For < u < 1 this metric gives S 2 the geometry of an oblate sphere, 
squashed along the direction of ipoo. For oj > 1, the metric is singular, changing from Rie- 
mannian to Lorentzian in a strip around the equator (orthogonal to ipoc)- Consequently, the 
pseudoenergy is no longer bounded below for u > U\ — 1 which, as we will see, has strong 
phenomenological consequences. 

The third term of F w is just the usual Faddeev-Skyrme term (quartic in spatial derivatives). 
The fourth and fifth terms together can be interpreted as a deformed potential 

U u {if>) = U(ip) - ia; 2 |^oo x ip\ 2 . (2.9) 

Since U attains its minimum at ip^ and is rotationally invariant, its hessian about ipoo must 
be /x 2 (-, ■) for some constant \i > 0, interpreted physically as the mass of mesons in the field 
theory. Hence, if cu > u) 2 = /U, F u is again unbounded below. A particularly convenient choice 
for U is 

U(ip) = ^ 2 (l-ip 2 3 ) 2 . (2.10) 
Then = (0, 0, 1) and the deformed potential is 

U UJ {ip)= l -{pi 2 -uj 2 ){l-iplf. (2.11) 
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This is the potential we use in all our numerical simulations. 

Since the model is invariant under global rotations of (p about ipoo, it has an associated 
conserved Noether charge, called isospin 

J= f {f (^x^) + (d(^- V ),^(^))}. (2.12) 

JM 

For uniformly isorotating fields of the form (2.5), this equals 

J = A(ip)uj (2.13) 

where the moment of inertia is 

A(V0=/ {hAocX./f + ld^-^l 2 }. (2.14) 

JM 

Hence F w {ijj) = V($) - \uo 2 k(ip) } while the total energy of the field (2.5) is V + T = V(i/j) + 
J 2 / (2A(ip)) . There are two natural variational problems for ip: 

1. For fixed u, extremize F u (ip) = V(ip) — \uo 2 k{ip); 

2. For fixed J, extremize Ejty) = V(i/>) + J 2 / (2A{ip)) . 

It is clear that these two problems are precisely equivalent: if ip solves 1, then it solves 2 with 
J = A(ip)uo, and if ip' solves 2, it solves 1 with to = J/A(i/j'). Previous studies [3] of isorotating 1 
solitons have used formulation 2, whereas in this paper we will mainly use formulation 1. This 
has several advantages. First, it is directly clear, as we have shown (using the Principle of 
Symmetric Criticality), that solutions of 1 correspond via (2.5) to genuine solutions of the 
field theory. Second, the Euler-Lagrange equation corresponding to 1 is a PDE, similar in 
structure to the static field equation of the Faddeev-Skyrme model, whereas the equation 
corresponding to 2 is a rather more complicated differential-integral equation. Consequently, 
it is a fairly simple matter to adapt existing numerical techniques, developed for the static FS 
model, to deal with problem 1. Third, formulation 1 makes it clear that, in the case fj, > 1, 
there is no reason why isospinning solitons should persist for frequencies uo G since 
is unbounded below when uo > min{l,/i}. Hence, we have the possibility that isospinning 
hopfions are destabilized by nonlinear velocity terms in the field equation before they reach 
the upper limit to — fj,. 

To conclude this section, we note that a general time dependent solution of the model 
tp(t, x) = R(ut)ijj(t, x) conserves both total energy T + V and isospin J, and hence conserves 

T + V - to J = F u {ij>) + lf (H 2 + \ H n\ 2 ). (2.15) 

Hence a local minimum of F u (xfj) corresponds to an orbitally stable solution of the model, 
since a small perturbation of ip is trapped close to ?/> by conservation of T + V — u J. So the 
solutions found by our numerical method, which can only find minima of F u , are guaranteed 
to be orbitally stable. One should note that conservation of T + V — to J does not directly 
imply that saddle points of F u are unstable. 

Actually, [3] concerns spatially rotating Skyrmions, but within the axially symmetric ansatz used therein, 
rotation is equivalent to isorotation. 
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3 An elastic rod model 



Before embarking on numerical investigation of the variational problems posed in the previous 
section, it is useful to consider analogous problems in a simple effective model of Faddeev- 
Skyrme solitons. A model based on elastic rods was introduced in [7], which successfully 
captures the qualitative features of static solitons with small Hopf degree on both M = M 3 [7] 
and M = M? x S 1 [6]. Although originally formulated only in flat 3-space, the model has a 
geometrically natural extension to any curved Lorentzian 4-manifold representing space-time, 
allowing one to model time- dependent solitons. 

The degrees of freedom in the elastic rod model consist of a smooth map x from an oriented 
2-manifold T into a 4-manifold Ai equipped with a Lorentzian metric g, and a section m of the 
normal bundle NT of unit length. The image of x is taken to represent the preimage under ip 
of the point — ipoo in S 2 antipodal to ipoo, and the section m represents the projection onto NT 
of the pull-back under if of a fixed vector in the tangent space at — ^oo. The SO (2) internal 
symmetry group of the Faddeev-Skyrme model, which rotates the target 2-sphere leaving ifi^ 
fixed, can be identified with the SO(2) structure group of NT,. The latter acts naturally on 
m, since the sphere bundle and the bundle of linear frames for NT can be identified. 

The elastic rod model is defined by an action comprising three terms which describe the 
effects of stretching, bending, and twisting of the rod. The stretching term S$ is simply the 
area of E, as measured using the metric 7 on T induced from g: 

S s (x) :=- [ Vol,. (3.1) 

This is nothing other than the Nambu-Goto action of string theory. 

The bending term Sb is constructed from the second fundamental form I of T. The latter 
is a section of Sym (T*T) ® NT and is defined by the equation, 

I(X, Y) = n N (V^°Y) VX, Y G T(TT). (3.2) 

Here V LC denotes the pull-back of the Levi-Civita connexion on M. to x*(TM) and tin 
denotes the orthogonal projection from x*(TA4) to NT. Since there are two naturally defined 
quadratic forms on Sym 2 (T*S), there are two ways to build an action from I. Introducing 
coordinates a a on T and on Ai, the two possible actions built from I = I^ b da a da b d^ are: 

Sb(x) := [ VTVttH and S' B (x) := [ g^-fXRVoL,. (3.3) 

These correspond to taking either the square of the trace or the trace of the square of the 
second fundamental form. 

While in general one should consider actions involving both terms Sb and S' B , in appli- 
cations in which M. is flat one may work solely with Sb without any loss of generality. The 
reason is that in this situation the equations of motion for both terms agree, as we now prove. 
Suppose that the Riemann curvature tensor for Ai vanishes. Then the Riemannian curvature 
R of 7 may be expressed in terms of the second fundamental form using Gauss's formula: 

iW = M«-«J- (3-4) 
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Contracting with ^ ac ^ hd and integrating over E results in the equation, 

J ScVol 1 = S B - S' b , (3.5) 

in which Sc denotes the scalar curvature of 7. The left hand side of this equation is equal 
to a linear combination of the Euler characteristic of E and an integral over the boundary 
of E, by the Gauss-Bonnet theorem for manifolds with indefinite signature [9]. In particular, 
it is stable to small local perturbations of x and thus has a trivial Euler-Lagrange equation. 
Therefore the Euler-Lagrange equations for Sb and S' B coincide. 

The twisting term St is built from the derivatives of m. Recall that a connexion V N can 
be defined on the normal bundle XE by orthogonally projecting the Levi-Civita connexion on 
x*TM: 

V^m := 7T N Vx°m VX eTJ:,me XE. (3.6) 

The twist rate 9 associated to the unit normal vector m is the projection of the covariant 
derivative of m to the subspace of XE orthogonal to m. More precisely, 9 is a section of T*E 
defined by the equation 

9(X)Vol 7 (Y : Z) = Vol g (V x m, m, Y, Z) VX, Y, Z e T(TE). (3.7) 

The twisting term is then 

Sr(x,m):= [ ^ ab 9 a 9 h Vol 1 = [ 9 A *~ f 9. (3.8) 
Jt, Jt, 

Equivalently, the twisting term may defined to be the L 2 norm of V N m. 
The complete action for the elastic rod model takes the form 

S(x, m) = AS s (x) + BSb(x) + CS T (x, m), (3.9) 

with A, B, C positive real parameters. Of these, two correspond to choices of units while only 
the dimensionless ratio B/C is non-trivial. 

3.1 Isospinning rods 

From now on the target manifold M. is taken to be M 3 x equipped with its standard 
Lorentzian metric, and E is taken to be compact. Let Y u denote the set of C 2 pairs (x, m) 
consisting of an immersion x : E — y M. and a section m of XE C x~ l TA4. For each m denote 
by m' the unique section of XE such that m, m' is a positively oriented orthonormal basis, 
and let K[ a ] denote the action of translation along S 1 ^ by a/u on R 3 x S^. The action of S 1 
on Y u corresponding to the action on X u is given by 

(x, m, a) t-j- (x[ a ], m[ a ]), X[ a ] :— K [ a ] o x, m[ a i := mcosa + m'sina. (3.10) 

The set of fixed points of this action will be denoted Y*j , and elements of this set may be 
described as follows. Let T be a compact 1-manifold, let x : T — )■ M 3 , and let m(s) be a unit 
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vector perpendicular to <9 s x. Suppose that a local coordinate s on T has been chosen so that 
\d s x.\ = 1; then the map is said to be arclength parametrised. Any element of Yj takes the 
form 

E = S^xT (3.11) 
x(t,s) = (t,x(s)) (3.12) 
m(t,s) = coswtm(s) + sinu;it(s) x m(s) (3.13) 

for some T, x, m. 

Recall that in the standard Frenet formalism for curves the unit tangent vector is defined 
to be t(s) := 9 a x(s), and the curvature is «(s) := \d s t\. The action for an element of Y*j is 
S = —{2n /uj)F u1 where 

F w (T,x,m):= ! (A - Ccu 2 + Bk 2 + C(t ■ <9 s m x m) 2 ) ds. (3.14) 

The critical points of the restriction of F u to Y^f 1 are critical points of S within Y u by the 
Principle of Symmetric Criticality [13]. 

As in the Faddeev-Skyrme model, the pseudoenergy for the rod model is bounded from 
below only for a range of values of u. Thus F u can only reasonably be expected to have critical 
points if u 2 < A/C, and attention will be restricted to this range. 

A crucial property of F u , which underpins the subsequent analysis, is that F w is equivalent 
to F up to a rescaling of lengths and energies: 

F w (T,x,m) = y/l-u 2 C/AF h, y/l - u 2 C/Ax, m) . (3.15) 

This means that properties of F u and its critical points can be gleaned directly from properties 
of Fo with no extra effort, in contrast with the full Faddeev-Skyrme model. 

Within the elastic rod model there is a conserved energy E and a conserved isospin J = 
—2C J r *8 which are the Noether charges associated with the symmetries of time translation 
and rotation of m. For elements of l^f 1 these take the form Ej = V + J 2 /2A and J = uA, 
with 

V = I (A + Bk 2 + C(t ■ d s m x m) 2 ) ds and (3.16) 

A = [ 2Cds. (3.17) 
Jr 

In particular, V is the energy functional for rods introduced in [7] (to see this, note that if 
m = sina(s)n + cosa(s)b then t • d s m x m = d s a — r with r denoting the torsion). 

Previous studies [7, 6] have found that the static rod model defined by V provides a useful 
approximation to Skyrme-Faddeev solitons if a value of approximately 0.85 is chosen for the 
dimensionless ratio C/B. The same choice should be equally useful in the Faddeev-Skyrme 
model with U ^ 0, since qualitative properties of solitons seem unaffected by the introduction 
of U [5]. In practice it is necessary to impose constraints on elastic rods in order to prevent 
self-intersections, but these can safely be ignored for the purposes of the present discussion. 
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The static rod model successfully captures qualitative features of all low-charge Faddeev- 
Skyrme solitons, with a minor exception at charge 4. The problem at charge 4 is that the 
minimal energy soliton is an .4.2,2 knot consisting of two superposed charge 2 solitons, and 
the rod model is unable to describe such a configuration. The charge 4 .4.2,2 soliton is the 
only known energy-minimising soliton which is a superposition of lower charge solitons, so 
this shortcoming of the rod model is very mild. Thus we expect that useful insight into the 
variational problems 1 and 2 posed in section 2 can be inferred by studying the corresponding 
problems for rods. 



3.2 Infinite strings 

Before considering solitons on M. 3 within the rod model, we will first discuss a simpler situtation 
of solitons with cylindrical symmetry. Let P > and let us replace physical space M = M 3 
by M = R 2 x S] n/P . The ansatz, 

T = Sl n/P , x(s) = (0, 0, s), m(s) = (cos(27rs/P), sin(27rs/P), 0), (3.18) 

describes the unique fixed point of a certain action of S 1 x SO(2) on . It is a critical point 
of the pseudoenergy at which the potential energy and moment of inertia take the values 

V = AP+ ^ ° and (3.19) 
A = 2PC. (3.20) 

Consider an extension of variational problem 2, in which Ej is extremized with respect to 
variations in P, x and m while J is held fixed. For each J the expression Ej = V + J 2 /2A has 
a unique minimum amongst S 1 x SO(2)-symmetric configurations. The values E (J), Pq(J) 
and ujo(J) of E, P and u at the minimum are easily obtained from equations (3.19) and (3.20) 
and the expression u = J/A: 



(4nyAC + P* (3.2i; 



P (J) = 27TW- + - — - (3.22) 




u (J) = ■ (3.23) 

A 

C 



{Att) 2 AC + P 2 



These equations show that both Eq(J) and Po(J) grow linearly with J for large values of J, 
and that u is approximately constant for large values of J. They also imply that the dimen- 
sionless quantity EoP^duo/dJ evaluates to in 2 at J = 0. These qualitative and quantitative 
predictions will now be tested in the Faddeev-Skyrme model. 

Imposing S 1 x SO (2) symmetry in the Faddeev-Skyrme model results in the following 
ansatz for ip in terms of a profile function f(r): 

i> = (sin f(r) cos(# + 2ttz/P), sin/(r) sin(# + 2nz/P), cos f(r)) . (3.24) 
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The expressions for V and A are 



V 



A 



7T 



(/ 



'\2 




i>oo 

2tt / (1 + (/' ) 2 ) sm 2 /rdr. 
Jo 



sin 2 / (1 + {f') 2 ) + [i 2 sin 2 / \ rdr (3.25) 



(3.26) 



The boundary conditions /(0) = and /(r) — > ir as r — > oo are imposed in order to maintain 
continuity at the origin and finiteness of V and A. 

Variational problem 2 may now be studied using the following numerical scheme. Let 
V(u),P) and A(u,P) denote the values of V and A at a critical point of = V — \oo 2 A. 
The values of the functions V(oo, P) and A (a;, P) at any pair (oo, P) may be determined by 
first determining a numerical solution / to the Euler-Lagrange equation associated with F w 
and then numerically evaluating the integral expressions for V and A. A minimum E (J) 
of E = V + |w 2 A with J = Au held fixed can be determined by numerically solving the 
J-preserving gradient flow equation for E in u and P: 



p = -d P E{u, P) + Md P J(u, P) 
Co = -d w E(u, P) + McU(w, P) 
dpEdpJ + d u Ed u J 



M 



dp J dp J + d^Jd^J 



(3.27) 





Figure 1: Infinite strings in the Faddeev-Skyrme model solving variational problem 2. The 
graphs show the energy E and period P as functions of the isospin, and also the energy E as 
a function of the parameter oo. 

The results of this procedure are presented in figure 1. The numerical results confirm the 
three qualitative predictions derived in the rod model for large J, namely that Eq(J) and 
Po{J) grow linearly with J, while ooq(J) tends to a constant. From the numerical data we find 
that 

2 du; 



EoPt 



dJ 



0.985 x 4tH 



(3.28) 



j=o 



This is within 2% of the value predicted using the rod model. 

To compare the rod and Skyrme-Faddeev models over a range of values of J we must 
choose values for A and C. We will take 



(4tt) 2 AC= 1.1, 



A 

c = 3 ' 5 ' 



(3.29) 
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The quantity (4ir) 2 AC has the dimension of energy squared, and its value has been chosen to 
fit the data for the infinite string. The quantity A/C has the dimension of length to the power 
of minus two. Its value has been chosen to give a reasonable fit for both the infinite string 
and for the solitons discussed below . The functions E Q (J), Po(J) and uj (J) in the rod model 
have been plotted in figure 1 alongside those for the Skyrme-Faddeev model. It can be seen 
that the curves for E (J) are in excellent agreement; the curves for uq and Pq are reasonably 
close for small J, but less close for large values of J. 

A final remark concerns the limiting value of uj. Previously we showed that the Faddeev- 
Skyrme pseudo-energy is unbounded below for u > 1, suggesting that the pseudo-energy 
does not have critical points for ui > 1. Our numerical results for infinite strings confirm 
that the range of u is bounded, however, the maximum value is around 1.4, and exceeds the 
expected value of 1. Thus for 1 < uj < 1.4 the pseudo-energy has critical points, despite being 
unbounded from below. 



3.3 Solitons 

Now let us consider the effect of isospin on elastic rods in M 3 . Suppose that (x , m ) is a 
critical point of ^=0 = V and that ^(xo, mo) = Vq. The moment of inertia Ao for this critical 
point may be determined by means of a virial theorem: the energy of a rescaled configuration 
(Ax , m ) is 

V(Ax , mo) = A / Ads + \ (Bk 2 + C(t • m x d s m ) 2 ) ds . (3.30) 
Jy A Vr 

Since (x ,m ) is a critical point the derivative of the left hand side with respect to A must 
evaluate to zero at A = 1, so the two integrals on the right hand side of the equation must be 
equal. It follows that 

A = jV . (3.31) 
Since F^ is related to F by equation (3.15) there exists a corresponding critical point 



(x /y 1 — cu 2 C / A,m ) of F^, and indeed all critical points of F^ are so obtained. The value 
of Fu at the critical point may be deduced from equation (3.15) to be 

F^ = y/1 - uj 2 C/AV . (3.32) 

Since A is proportional to the length of x(T), the value of A at this critical point is 

A = - 1 ^V . (3.33) 

It follows that J = uj A and E = F^ + u 2 A take the values 

J = ; UC/A V and (3.34) 

y/1 ~ UJ 2 C/A 

E = Vp (3.35) 

V 2 + J 2 A/C (3.36) 
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at the critical point. 

The preceding formulae contain useful qualitative information about the solutions to vari- 
ational problems 1 and 2 in the elastic rod model. Equation (3.32) implies that the values of 

at its critical points all scale with u in the same way. In particular, the global minimum at 
oj = remains the global minimum of F u for all other values of u. As u approaches a critical 
value the moment of inertia, angular momentum, energy and size diverge. 

The behaviour of solutions to the second variational problem is similarly uniform: equation 
(3.36) captures the dependence of the critical values of Ej on J. The global minimum at 
J = remains the global minimum for all values of J, because the right hand side of (3.36) 
is monotonic in Vq. At large values of J the energy, moment of inertia and size grow linearly 
with J, while u tends to a finite constant. Similar behaviour should be expected from solitons 
in the Faddeev-Skyrme model. 



4 Numerical results 

In this section the results of numerical simulations of the Faddeev-Skyrme model will be 
presented. Local minima of the pseudo-energy functional F w (xjj) have been found for a range 
of values of the Hopf degree Q and of u. The majority of simulations were carried out with 
« = 2. This choice guarantees that oj\ < 0J2, and thus allows us to investigate the possibility 
that solitons become unstable to processes other than pion decay. Other values of \x were also 
investigated with Hopf degrees one, two and three. 

The numerical algorithm employed was similar to that used in [8]. A simple forward 
differencing scheme was used on a cubic lattice with a lattice spacing Ax = 0.1 on a grid 
typically containing (240) 3 latice points. Well-chosen initial configurations were evolved using 
the quasi-Newton BFGS method to minimize the pseudoenergy functional F^i^i). Simulations 
were considered to have converged to local minima if the sup norm of the gradient of F u (i/j) 
was less than 0.01 Ax 3 . 

Each of our simulations began at u = and proceeded by making small increments in u. 
At each value of u the algorithm was allowed to converge to a local minimum of the pseudo- 
energy, and this local minimum was then used as the initial configuration for the next value 
of u. Initial configurations at u = were created using Sutcliffe's rational map ansatz [15]. 
This involves identifying physical space M 3 with the three-sphere S 3 C C 2 using a map 

(Z l5 Z ) = ((*i + ^x 2 ) S ^l, cos f(r) + ix a ^Q) , (4.1) 



in which f(r) is some monotonic function of r = \/ x\ + x\ + x\ satisfying the boundary 
conditions /(0) = 7r, /(oo) = 0. Initial configurations are given by maps of the form 



V>i + #2 z?z ( 



W = = (4-2) 



1 + ?/>3 Zf + Z, 



with a,b,a,(3 natural numbers. These maps have Hopf degree Q = ab + /3a, so that an 
initial configuration with any desired Hopf degree can easily be constructed. The largest value 
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of ui that we were able to explore was to = 1, reflecting the fact that the pseudo-energy is 
unbounded below thereafter. However, the existence of saddle points of F u with uj > 1 is not 
ruled out. 

Below we present location curves and energy density isosurfaces for solitons, and graphs of 
their energies, both as functions of uj and as functions of isospin. As in recent studies [15, 5], 
the location curves plotted are the tubular preimages of the curve ^3 = —0.90. 

Energy curves obtained from the rod model are plotted alongside those obtained from the 
numerics, in order to ascertain whether soliton energies follow the expected behaviour. The 
rod model curves are obtained from equations (3.35) and (3.36). The value of A/C is taken to 
be 3.5 for all solitons with \i — 2, and for each soliton the value of Eq is taken to be the soliton 
energy at uj = 0. In graphs containing energy curves for several solitons, we choose to plot 
only one or two curves from the rod model, in order to avoid cluttering the graph. Graphs are 
presented using normalised units of energy in which E = E/(16ir 2 V2) and J = J/(16ttV2); 
with this convention Ward's conjectured lower bound on the energy [16] is E > Q 3 ^ 4 . 

In the previous section it was argued using the elastic rod approximation that the sizes and 
energies of solitons should scale in a uniform way both with u and with the isospin J. These 
predictions are expected to hold for relatively small values of uj or J, but at larger values some 
deviation from this approximation should be expected. It is well-known that for most values 
of Q the Faddeev-Skyrme energy has several local minima with similar energies. Thus it is 
plausible that at large values of J or uj small deviations from the elastic rod approximation can 
occasionally result in the crossing of energy curves, and hence in a change in the qualitative 
features of the global minimizer of energy. The results presented below confirm that the E(u) 
curves follow the predictions of the rod model at small uj. Somewhat surprisingly, the E(J) 
curves continue to do so even at large J. We also observe occasional crossings in both the 
E(J) and the £"(0;) curves. 



4.1 Degrees one, two and three 




Figure 2: Energy (in units of Q 3//4 ) as a function of angular frequency uj and as a function of 
J for solitons with Q = 1,2,3 and fi = 2, and energy curves derived within the rod 



isospin 
model. 
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lAx,x 2A 2 ,i 3A\i 

Figure 3: Position curves for solitons with charges \ <Q<3&tfi = l. The light green and 
blue position curves correspond to u = and u — 1, respectively. 




0.2 0.4 0.6 0.8 1 1.2 0.5 1 1.5 2 2.5 3 3.5 

U) J 

Figure 4: Energy as a function of angular frequency u and as a function of isospin J for 
Q = 1, 2, 3 for a range of values of fi. The transition 2.42,1 — > -4.2, i at /i = 1 and u = 0.86, is 
labelled. 



For Hopf degrees 1, 2 and 3 simulations were carried out with a range of values of \i. The 
energy curves for /i = 2 are plotted in figure 2 alongside the predictions from the rod model. 
These clearly show that the soliton energies scale as expected. 

In figure 4 soliton energy curves for a range of values of \i are plotted. Figure 3 shows the 
soliton location curves for \i — 1 at a; = and u — 1. For most values of u and fi the preferred 
configurations are of type Ai t i, 2A 2 ,i and 3.43,1, as in the static case. Here, and in the sequel, 
we are using the notation introduced by Sutcliffe to label Hopf soliton shapes [15]. Briefly, 
QA n ,m denotes an axially symmetric hopfion of charge Q, where the position curve (preimage 
of — ipoo) is a circle, and the preimage of a regular value close to — ipoo is a disjoint union of 
m closed curves, each winding n times around the circle. In fact Q = nm, so including Q 
in the label is redundant, but convenient. A soliton of type QA n ^ m has the same qualitative 
form, but with axial symmetry weakly broken, so the position curve is not exactly circular. 
Later we will encounter solitons whose position curves are links of two components. These 
will be denoted QCJ^\ where the subscripts denote the Hopf charges of each component, and 
the superscripts denote the extra Hopf degree of each component due to its linking with the 
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other component. Again, we include Q, though this is redundant (since Q = p + q + a + b). 
We will also encounter hopfions whose position curves are torus knots of type (a, b) (where a 
and b denote the windings of the curve around the S 1 factors in T 2 ). We denote these Q/C a ,£> 

When Q — 2, \i — 1 and uj > 0.86, our algorithm converges to a buckled 2^4 2 ,i rather than 
an axial 2^4 21 configuration. This configuration did not appear in the earlier study [1], whose 
attention was restricted to axial symmetry. The buckled configuration seems to be a local 
minimum of the pseudo-energy, so should be orbitally stable (as explained in section 2). It is 
likely that an axial 2.4.2,1 configuration continues to exist above u = 0.86, albeit as a saddle 
point rather than a minimum of F^. Since our algorithm only finds local minima and not 
saddle points we are unable to determine whether the 2^4.2,1 has a greater or smaller energy 
than the 2.4.2,1 for any given value of u or J. 

When fi < 1 our simulations terminate at u = fi, as expected. When fi > 1 our algorithm 
ceases to find any critical points when uj = 1 (although in a few cases it is possible to go 
slightly beyond oj — 1). This is again consistent with our expectations, as the pseudo-energy 
is not bounded from below when oj > 1. However, the existence of solutions with oj > 1 is not 
ruled out, as they may continue to exist as saddle points of the pseudo-energy. 

The borderline case fi = 1 is particularly interesting: the graphs of E(uj) grow rapidly as 
oj approaches 1. The rod model predicted that soliton energies would diverge at some finite 
value of oj, and our numerical results suggest that for /i = 1 this value is very close to 1. 

4.2 Degree four 



44 4 ,i 4A, 2 /4£!;1 saddle 

Figure 5: Position curves for solitons with Q = 4 and fi = 2. In the left and middle pictures 
the light green and blue position curves correspond to u = and uj = 1 respectively. The 
right picture is an intermediate state between 4^4 2 ,2 and ^C\\ at u = 0.58. 

At degree 4 we have been able to find configurations of types 4*4 2 ,2, 4*4 4j i and 4:C\'\; these 
are depicted in figure 5. The 4*4 2 ,2 is axially symmetric and may be thought of as two adjacent 
2.4 2 ,i solitons. The position curve of the 4^4. 2 ,2 in our model consists of two adjacent circles, 
whereas in the massless Faddeev-Skyrme model it is a single circle. 

The graphs of energy as a function of u; or J are shown in figure 6. The 4A^i configuration 
exists for all values of u in the range [0,1], but the 4^4 2 ,2 and 4:C\'\ configurations could only 
be found in the ranges [0,0.65] and [0.57,1] respectively. At u = 0.57 the 4«4. 2 ,2 and 4:C\'\ 
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Figure 6: Energy as function of the parameter oj and as a function of isospin J for solitons 
with Q = 4 and \i — 2. In the right panel we only present data in the range 0.4 < J < 1.1 to 
better illustrate the energy crossings. 



configurations are degenerate in both energy and pseudo-energy, from which we infer that 
the 4^4.2,2 critical point undergoes a bifurcation at this point. We were also able to find a 
configuration which seems to be a saddle point of the pseudo-energy at oj = 0.58, suggesting 
that the bifurcation is a pitchfork bifurcation. This saddle is illustrated in figure 5 but its 
energy has not been plotted. 

When oj > 0.65 our algorithm is unable to find the 4.4.2,2 and instead converges to the 
4£^'}. Since the 4^2,2 configuration is axially symmetric it is likely to continue to exist as 
a critical point of the pseudo-energy beyond oj = 0.65; the most likely reason for our code's 
failure to find it is that is a saddle point rather than a local minimum. 

Both the E(oj) and the E(J) graphs exhibit crossings. The minimum of the energy at 
J = is the 4.4.2,2 configuration. All three energy graphs cross in the range 0.5 < J < 1.0, 
and for J > 1.0 the energy-minimzer is the 4:C\'\ link. It is likely that the energy-minimizer 
for fixed oj is always the 4^2,2, but for oj > 0.65 it is difficult to be sure of this without being 
able to construct the 4^4 2 ,2 configuration. We note that, despite the presence of crossings, the 
energy curves generally follow the rod model predictions, the only significant deviant being 
the E(oj) curve for the link when oj > 0.6. 

4.3 Degrees five to eight 

At degree five we found two distinct local minima of the pseudo-energy, namely a 5C\'1 link 

and a 5^4.5^ buckled ring, see figure 7. When oj = the link has the lower energy. The two 
curves for E(oj) cross when oj = 0.82 and thereafter the ring has the lower energy. In contrast, 
the E(J) curves for these configurations do not cross and the link has the lower energy for 
any given value of J. 

At degree six we again found two distinct local minima of the pseudo-energy. These were 
links of type GCl'l and 6^3'} , shown in figure 8. Unlike the study [5], we did not find a 6A^2 
configuration; this could be due to a different choice of potential function or a different choice 
of the mass parameter [i. The ^C^X nn k has a lower energy than the 6/^3 \ link for all values 
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5*4-5,1 5jC^'2 




Figure 7: Solitons with Q = 5 and /i = 2. Top row: position curves for the 5^5^ and 5£\'l 
configurations, with light green and blue position curves corresponding to u = and uj = 1 
respectively. Bottom row: energy as function of the parameter uj and as a function of isospin 
J. 

of uj and for all values of J. 

At degree seven the only minimum of the pseudo-energy found was a 7/03,2 knot, and this 
is shown in figure 9. 

At degree eight the three energy minima found were a link of type 8C^'\ and a knot of type 
8/C 3i2 , see figure 10. We did not find a soliton corresponding to the 8^4,2 configuration in [5]. 
Within the limits of numerical accuracy the 8C^' 3 and 8^2 configurations are degenerate in 
energy when uj = 0. Asw increases the knot energy grows faster than that of the link, so that 
the link has the lower energy. When uj reaches the value 0.38 the link collapses to the knot, 
which has a smaller pseudo-energy. Our algorithm is unable to find the link at larger values of 
uj. It is likely that when uj > 0.38 the link is an unstable critical point of the pseudo-energy, 
and that it continues to have a lower energy than the knot. As far as we can tell, the E(J) 
curves for the two solitons coincide, at least within the range of values of J in which both can 
be constructed. Note that, since their energies are indistinguishable at ou — 0, this is precisely 
what the rod model predicts. 
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Figure 8: Solitons with Q = 6 and \x = 2. Top row: position curves for the QC^2 and 6jC 3 '\ 
configurations, with light green and blue position curves corresponding to u = and uj = 1 
respectively. Bottom row: energy as function of the parameter uj and as a function of isospin 
J. 

5 Conclusions 

To summarize, we have reformulated the problem of constructing internally spinning Hopf 
solitons of frequency u as a variational problem for pseudoenergy F u , and solved this problem 
numerically for Hopf charges 1-8, < uj < 1 without imposing any spatial symmetries. 
For each local minimizer of the static energy F we have constructed a curve of isospinning 
hopfions, parametrized by uj, for which the total energy E{uj) and conserved isospin J(uj) 
were calculated. Generically, the solitons persist for all uj < min{l,/i}, where /i is the meson 
mass of the model, and their qualitative shape is independent of uj. We noted two exceptions 
to this: for Q — 2, fi = 1, the solitons undergo a bifurcation as uj increases through 0.86, 
branching from 2*4.2,1 to 2.4.2,1 (he. they lose axial symmetry); for Q = 4, fi = 2, the axially 
symmetric 4^4. 2 ,2 soliton branches to a link 4£}'} as uj increases through 0.58. Even when no 
such bifurcation occurs, the shape of the lowest energy soliton with fixed uj can change as uj 
varies, because the graphs of E{uj) for different solitons with the same Hopf degree can cross. 
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UJ J 

Figure 9: Solitons with Q = 7 and \i — 2. The light green and blue curves in the left figure 
are position curves for solitons at u = and ou — 1 respectively. The two graphs show energy 
as a function of the parameter u and as a function of isospin J. 




We found that this happens for Q = 5, [i — 2, for example. Similarly, the shape of the lowest 
energy soliton with fixed isospin J can change as J varies, if the graphs of E(J) for different 
soliton branches cross, as happens, for example for Q = 4, \i — 2. Such crossings seem to be 
comparatively rare, however. 

We have developed a time-dependent extension of the the elastic rod model of Hopf solitons, 
which predicts simple, uniform scaling behaviour of the total energy and isospin of hopfions 
as bj varies, namely 

h{(jj) = J{u}) = = (5.1) 

where Eq is the static soliton energy and is an (unknown) parameter, depending only on 
the details of the Faddeev-Skyrme model under consideration (but independent of Q, and of 
which branch of local minimizers one is following). Despite having only one free parameter, 
the rod model fits the E(oj) graphs well for small u, and the E(J) graphs well for all J. The 
rod model predicts that the energy graphs never cross, and that solitons never change shape 
(they just dilate uniformly), as to (or J) varies which, again, is in good qualitative agreement 
with our numerical results. 
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Figure 10: Solitons with Q = 8 and /x = 2. The top row shows position curves at w = (light 
green) and w = 1 (blue). The bottom row shows energy as a function of the parameter u and 
as a function of isospin J. 
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